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ABSTRACT 

The AKARI FU-HYU mission program carried out mid-infrared imaging of several well studied Spitzer fields preferentially selecting 
fields already rich in multi-wavelength data from radio to X-ray wavelengths filling in the wavelength desert between the Spitzer IRAC 
and MIPS bands. We present the initial results for the FU-HYU survey in the GOODS-N field. We utilize the supreme multiwavelength 
coverage in the GOODS-N field to produce a multiwavelength catalogue from infrared to ultraviolet wavelengths, containing more 
than 4393 sources, including photometric redshifts. Using the FU-HYU catalogue we present colour-colour diagrams that map the 
passage of PAH features through our observation bands. We find that the longer mid-infrared bands from AKARI (IRC-L18W 18 
micron band) and Spitzer (MIPS24 24 micron band) provide an accurate measure of the total MIR emission of the sources and 
therefore their probable total mid-infrared luminosity. We also find that colours incorporating the AKARI IRC-S 1111 micron band 
produce a bimodal distribution where an excess at 1 1 microns preferentially selects moderate redshift star-forming galaxies. These 
powerful colour-colour diagnostics are further used as tools to extract anomalous colour populations, in particular a population of 
Silicate Break galaxies from the GOODS-N field showing that dusty starbursts can be selected of specific redshift ranges (z=1.2 - 
1.6) by mid-infrared drop-out techniques. The FU-HYU catalogue will be made publically available to the astronomical community. 

Key words. Infrared: source counts, Surveys - Cosmology: source counts - Galaxies: evolution. 



1. Introduction 

Studies with the Infrared Space Observatory (ISO) of the 
Hubble deep fields, North and South (HDF-N, HDF-S) have 
revealed star formation rates at least comparable to or higher 
than those of the optic al/UV studies (Ma nn et al. (2002)| > . At 
submillimetre wavelengths, surveys in the same HDF-N field 
with SCUBA on the JCMT has revealed a large (>3000 deg" 2 
at S 850 >2 mJy) population of strongly evolving sources 
with bolometric luminosities > 10 12 Lq and star formation 
rates of -SOO-lOO OMoyr' 1 with a median redshift of 2.4 
(Chapman et al.(2005)[ ). The overwhelming conclusion is that 
the star formation rate at z~l-2 requires significant evolution 
in the IR galaxy population from the current epoch. Deep ob- 
servations with the Spitzer Space Telescope have confirmed this 
strong evolution in the galaxy population out to the intermediate 
redshift range probed by ISO (z~0.3-l) and furthermore pro- 
vided insight into the higher redshift Universe in the so called 
redshift desert z~l-3. To connect the local and intermediate 
redshift Universe to the higher and high z Universe observed 
by Spitzer and SCUBA, comprehensive multiwavelength imag- 



ing and spectroscopy is required throughout the extragalactic 
population. The Fields surveyed by the Spitzer space telescope 
are some of the most richest fields in multiwavelength data in 
the entire sky. The SWIRE wide area survey (50 square de- 
grees over 6 fields, Lonsdale et al.(2 004)[), the GOODS su rvey 
(300 square arcmins over 2 fields, |Giavalisco et al.(2004)) >, the 
other GTO fields and the Spitzer First Look Survey (FLS, e.g. 
|Frayer et al.(20 06)) as well as having coverage over the 3.5- 
160jum wavelength range also include a wealth of ancillary data 
at other wavelengths from radio to X-rays. However, Spitzer has 
limited imaging capability between 8yum< A <24pm and no 
spectroscopic capability shortward of 5pm leaving a conspicu- 
ous gap in these ranges in both wavelength and redshift. 

The AKARI satell ite, (formerly known as ASTRO-F, 
Murakami et al. (2007)[ ) is the first Japanese space mission ded- 
icated to infrared astrophysics and was launched on board 
JAXA's M-V Launch Vehicle No. 8 (M-V-8) at 6:28 a.m. on 
February 22, 2006 Japan Standard Time, JST (February 21st, 
9.28 p.m. UT) from the Uchinoura Space Center (USC) on the 
southern tip of Japan. The satellite is in a Sun-synchronous po- 
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Table 1. Photometric Filter band definitions for the three cam- 
eras of the IRC. 



Camera 


Array Area* 
(pixels) 


Pixel Scale 
(arcsec) 


Band 


(pm) 


(jim) 


NIR 


319x412 


1.46x1.46 


N2 


2.4 


1.9-2.8 








N3 


3.2 


2.7-3.8 








N4 


4.1 


3.6-5.3 


MIR-S 


233x256 


2.34x2.34 


S7 


7.0 


5.9-8.4 








S9W 


9.0 


6.7-11.6 








Sll 


11.0 


8.5-13.1 


MIR-L 


246x239 


2.51x2.39 


L15 


15.0 


12.6-19.4 








L18W 


18.0 


13.9-25.6 








L24 


24.0 


20.3-26.5 



'■ This is the effective imaging area of the array 



lar orbit at an altitude of 700 km and a period of 100 minutes. 
AKARI has a 68.5 cm cooled telescope with two focal plane 
instruments, namely the Far-Infrared Surveyor (FIS) and the 
Infrared Camera (IRC). The FIS has two 2-dimensional detector 
arrays and observes in four far-infrared bands between 50 and 
180 yum (Kawa daet al. (2007)| l. The IRC consists of three cam- 
eras covering 1.7- 26 pm in 9 bands (see Table [TJ with fields 
of view of approximately 10' x 10' ( |Onaka et al. (2007)| l. Both 
instruments have low- to moderate-resolution spectroscopic ca- 
pability. 

The AKARI mission is primarily a survey mission with an 
All-Sky Survey in the far-infrared and large mid-infrared legacy 
survey programs being carried out in the region of the North 
Ecliptic Pole (Matsuhara et al. (200 6)) and the Large Magellanic 
Cloud. However in addition to these Large Surveys (LS), a cam- 
paign of guaranteed time observations, referred to as Mission 
Programs (MP) was also undertaken. 

We report the initial results from the AKARI-FU-HYU 
Mission Program (Follow-Up Hayai-Yasui-Umai). The FU- 
HYU MP strategically targeted well-studied fields to maximise 
the legacy value of the AKARI data and in particular to target 
Spitzer fields in the above mentioned conspicuous wavelength 
"gap", at wavelengths of 11, 15 & 18jum. The FU-HYU obser- 
vations fill in vital gaps in the wavelength coverage (see Figure 
[TJ and provide invaluable insight into the connection between 
ISO & Spitzer populations and linking the far-infrared Universe 
to the high redshift sub-mm Universe as observed by the SCUBA 
instrument on the James Clerk Maxell Telescope. 

The FU-HYU program imaged three major fields (GOODS- 
N, Lockman Hole, ELAIS-N1, see Section 13 and in this work 
we report on the initial data reduction and results from the 
GOODS-N observations. Observation in the Lockman Hole will 
be reported by Serjeant et al. (2009) 

Unfortunately, the FU-HYU observations in the GOODS-N 
field were not dithered during the observation operation resulting 
in many hot and bad pixels leaving holes in the final co-added 
images. Moreover the distortion correction in the standard IRC 
pipeline toolkit blurs the hot pixels compounding the problem. 
In Section [3] we describe in detail our extensive additional pro- 
cessing utilizing intrinsic jitter in the spacecraft observations as 
an effective dithering and interrupting the standard IRC pipeline 
toolkit before the distortion correction stage to carry out inde- 
pendent re-gridding of the pixels onto a finer mesh, co-adding 
and re-binning, resulting in final dithered images without the hot 
pixels and line artifacts seen in the original processed data. 
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Fig. 1. The comprehensive coverage of the AKARI and Spitzer 
bands in the FU-HYU fields. The combination of the AKARI and 
Spitzer bands completely cover the features in the mid-infrared 
spectra of star-forming galaxies (STFG, shown in the figure at 
three different redshifts) locally, out to redshifts of ~2. 



The cross associations with other data sets in the GOOD-N 
field, most notably the Spitzer data and the construction of the 
FU-HYU catalogue is described in Section [4] The results and 
conclusions are given in Sections [5] & [^respectively. 



2. The AKARI FU-HYU program 

The objective of the FU-HYU Mission Program is to add AKARI 
data to contemporary data rich fields. On selection of our tar- 
get fields, the constraint on the visibility of the AKARI satel- 
lite had to be seriously considered. The orbit of AKARI is Sun- 
synchronous meaning that the ecliptic poles enjoy many point- 
ing opportunities whilst fields near the ecliptic may only be vis- 
ible on a few orbits every six months. Taking as an example the 
two well known fields of ISO-EL AIS Nl ( |01iver et al. (2000)| ) 
at an ecliptic latitude o f f3 -73° and the UKIDSS-UDS field 
( |Lawrence et al.(2007)| > at ft — 17°, we find that the former is 
observable of 17 orbits during half a year whilst the latter is ob- 
servable only 5 times. Given that not every AKARI orbit is avail- 
able for observation (calibration, directors time, etc) and further 
constraints placed on the observations by the necessity of choice 
of the deep imaging AOT (see below), there exist very clear cri- 
teria for our choice of target fields. 

For the FU-HYU program we selected three contemporary 
data rich fields; 

- The GOODS North field (GOODS-N) 

- The ELAIS-N1 field 

- The Lockman Hole 

The ELAIS-N1 field (R.A.= 16h09m20s; 
DEC=+54°57'00", J2000, /?=73°) is one of the the most 
richest fields in multiwavelength data on the entire sky. The 
ELAIS-N1 field was selected as part of the Extragalactic Large 
Area ISO Survey (ELAIS, piver et al. (2000)1 ). The fields were 
selected on the basis of their high ecliptic latitude (J3 >40°, to re- 
duce the impact of zodiacal dust emission), low cirrus emission 
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(hoOfim < 1-5 MJy/sr) and the absence of any bright (Sn^,,, > 
0.6 Jy) IRAS 12/im sources. A field of ~2 square degrees 
was observed to -1.0, 0.7, 70, 223mJy at 6.7, 15, 90, 175yum 
respectively and has extensive spectroscopic and imaging follow 
up at optical and near-IR wavelengths in the U, g, r, i , Z, J, H, 
K bands (Vaisanen et al. (2002)| |Gonzalez-Solares et al.(2005), 
Serjeant et al. (2004)| i. X -ray follow up with Chandra in the 



central region of Nl and a comparison with ROSAT has 
been made by |Manners et al. (2003)1 jWillott et al. (2003)| 
Basilakos et al. (2002) Radio observations also exist at 20cm 
qCiliegi et al.(1999)| l. The ELA IS Nl field was also partially 
observed with the Ha survey of Pascu aTet al. (2001)) Most sig- 
nificantly, ELAIS-N1 was a Spitzer SWIRE Field and was also 
selected as part of the GOODS Science Verification observations 
(Chary et al. (2004)). This is one of the deepest ever observa- 
tions of the 24pm sky down to around 20yuJy over 130arcmin 2 
centred on R.A.=16h09m00.16s, Dec=+54°55'50.9" and is the 
centre of the FU-HYU observation in the ELAIS-N1 field. Of 
the 3 FU-HYU fields, the ELAIS-N1 field has the best visibility 
and was observed with a total of 22 pointings over one 10'x 10' 
field of view (7 in the IRC S 1 1 band, 8 in the IRC L 1 5 band and 
7 in the IRC L18W band, see Table©. 

The Lockman Hole field (R.A.= 10h52m43s; 
DEC=+57°28'48", J2000, £=53°) has one of the lowest 
HI column densities and background cirrus emission on the sky 
(Iiou/im < 0.4 MJy/sr). The field has been covered by Spitzer 
by both IRAC and MIPS instruments ( |Le Floc'h et al. (2004)| 
Huang et al. (2004)). The Lockman field has also been mapped 



as part of the SCUBA Half Degree Extragalactic Survey 
(SHADES, |Mortier et al. (2005)| ) and recently supplemented 
by equivalently deep 1.1mm data with the AzTEC camera 
( |Austermann et al.(2009)| ). The Lockmann Hole is also the 
site of the UKIDSS DX S survey JHK data to K=21, H=22, 
J=22.5 over 8.75 sq deg. ( |Lawrence et al.(2007)). At radio and 
X-ray wavelengths there is VLA data at 1.4GHz to about 15yuJy 
(Big gs & Ivison(2006)| l and XMM data ( |Hasinger et al. (2001)| l 
respectively. The objective of the Lockman field observations 
was to gather AKARI data on the submillimetre sources in the 
field however the visibility was unfortunate in that the oppor- 
tunities to observe the Lockman Hole fell between observation 
phases in the AKARI mission planning (i.e. PV phase and the 
end of Phase I observations) limiting the number of possible 
10'x 10' fields of view to only three (see Table 0. The results 
from the FU-HYU Lockman Hole Observations are presented 
in |Serjeant et al. (2009)| 

The GOODS-N field (RA.= 12h36m 50s; 
DEC=+62°12'58", J2000, /3=57°, |Dickenson et al.(200T)) i has 
perhaps the best multi-wavelength coverage of any pencil-beam 
survey (160' 2 ), and is presented here as the first priority for our 
initial analysis. It is of course co-located with the Hubble deep 
field North ( fWiliiams etal. (1996)|l with optical data from the 
HST treasury programme (Giavali sco et al.(2004)[>, deep ground 
based UBVRIJHKz' Subaru/KPNO data (|Capak et al.(2004) i 
and Keck optical spectroscopy ( Reddy et al. (2006) I. The 



GOODS-N field has also been surveyed in the infrared with 
ISO ( [Serjeant et al. (1997)| l & Spitzer (See Section g2]l and at 
submillimetre wavelengths with the SCUBA and AzTEC in- 



struments at 850pm and 1 . 1mm respectively ( Pope et al. (2005) 



|Perera etal. (20 08)). The visibility of the GOODS-N field was 
marginal for AKARI and the emphasis was placed on a single 
10'x 10' field of view in the unique IRC L18W (3 pointings) 
and IRC N4, Sll (6 pointings each) bands. The full details of 
the observations are given in Table [2] Combining the AKARI 
and Spitzer data in the GOODS-N field together provides com- 
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Fig. 2. Coverage of the AKARI FU-HYU area within the 
GOODS-N field. The image is the Spitzer IRAC 3.6pm band 
image. The large round rectangle depicts the Spitzer coverage in 
the MIPS24yimi band. The larger of the two inner squares is the 
AKARI-IRC MIR-L L18W band coverage and the smaller inner 
square is the coverage of the IRC MIR-S SI 1 & NIR N4 bands 
(i.e. the foot print of the detector arrays) 



prehensive coverage of the entire near-mid infrared spectrum 
of dusty galaxies and their multitude of emission/absorption 
features from redshifts to 2. The supreme coverage of the 
combined AKARI and Spitzer bands is highlighted in FigureQ] 

In Figure [2] we show the FU-HYU imaging field of view 
(-10x10 arcmin for both the MIR-L and MIR-S/NIR IRC chan- 
nels) in the GOODS-N field. These, along with the Spitzer MIPS 
24-pm coverage are overlaid on the image in the Spitzer IRAC 
3.6yum band. 



3. Data reduction 

3.1. Standard processing with the IRC pipeline 

In nominal circumstances, AKARI-IRC data is reduced by 
the standard IRC imaging pipeline inside the IRAF environ- 
ment fLorente et al. (2007)[ ). The IRC pipeline consists of three 
stages; 

- Pre-pipeline (red-box): Takes the packaged raw data as in- 
put and then slices these into individual exposure frames and 
creates an exposure frame observation log. 

- Pipeline (green-box): The output of the green box produces 
the basic data for all exposure frames, processing each frame 
individually, correcting for instrumental effects (e.g. dark 
subtraction, cosmic ray rejection, saturation correction, flat 
fielding, distortion correction, etc). 

- Pipeline (blue-box): Co-add individual frames together and 
astrometry solution. The output is the final image and signal- 
to-noise map 

A series of three AOT sets for pointed observations with the 
IRC have been prepared. Two of these, referred to as IRC02 
and IRC03, are intended for medium deep multi-band surveys, 
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Table 2. Observation Log for the FU-HYU Mission Program 



Field RA 


DEC 


Observation ID 


Band 


AOT 


Parameters 


Date 


Status 


GOODS-N 189.2079 


62.2161 


1320028-001 


Sll 


IRC05 


c;N 


19/11/06 


Observed 






1320029-001 


Sll 


IRC05 


c;N 


19/11/06 


Observed 






1320030-001 


Sll 


IRC05 


c;N 


20/11/06 


Observed 






1320031-001 


Sll 


IRC05 


c;N 


21/11/06 


Observed 






1320032-001 


Sll 


IRC05 


c;N 


21/11/06 


Observed 






1320033-001 


Sll 


IRC05 


c;N 


21/11/06 


Observed 






1320025-001 


L18W 


IRC05 


a;L 


19/11/06 


Observed 






1320026-001 


L18W 


IRC05 


a;L 


19/11/06 


Observed 






1320027-001 


L18W 


IRC05 


a;L 


19/11/06 


Observed 


ELAIS-N1 242.2507 


54.9308 


1320226-001 


Sll 


IRC05 


c;N 


19/7/07 


Observed 






1320226-003 


Sll 


IRC05 


c;N 


19/7/07 


Observed 






1320226-004 


Sll 


IRC05 


c;N 


19/7/07 


Observed 






1320226-005 


Sll 


IRC05 


c;N 


19/7/07 


STT WDT error 






1320226-002 


Sll 


IRC05 


c;N 


19/7/07 


Data lost 






1320226-006 


Sll 


IRC05 


c;N 


20/7/07 


Observed 






1320226-007 


Sll 


IRC05 


c;N 


20/7/07 


Observed 






1320014-001 


L15 


IRC05 


b;L 


15/1/07 


Observed 






1320013-001 


L15 


IRC05 


b;L 


16/1/07 


Observed 






1320015-001 


L15 


IRC05 


b;L 


17/1/07 


Observed 






1320235-001 


L15 


IRC05 


b;L 


21/7/07 


Observed 






1320235-002 


L15 


IRC05 


b;L 


21/7/07 


Observed 






1320235-003 


L15 


IRC05 


b;L 


21/7/07 


Observed 






1320235-004 


L15 


IRC05 


b;L 


22/7/07 


Data lost 






1320235-005 


L15 


IRC05 


b;L 


22/7/07 


Observed 






1320232-001 


L18W 


IRC05 


a;L 


20/7/07 


Observed 






1320232-002 


L18W 


IRC05 


a;L 


20/7/07 


Observed 






1320232-003 


L18W 


IRC05 


a;L 


20/7/07 


Observed 






1320232-004 


L18W 


IRC05 


a;L 


21/7/07 


Observed 






1320232-005 


L18W 


IRC05 


a;L 


22/7/07 


Observed 






1320232-006 


L18W 


IRC05 


a;L 


22/7/07 


Observed 






1320232-007 


L18W 


IRC05 


a;L 


22/7/07 


Observed 


Lockman Hole 162.7016 


57.5555 


1320102-001 


L15 


IRC05 


b;L 


7/5/07 


Observed 






1320103-001 


L15 


IRC05 


b;L 


7/5/07 


Observed 






1320104-001 


L15 


IRC05 


b;L 


7/5/07 


Observed 


162.9408 


57.4731 


1320099-001 


L15 


IRC05 


b;L 


7/5/07 


Observed 






1320300-001 


L15 


IRC05 


b;L 


8/5/07 


Observed 






1320301-001 


L15 


IRC05 


b;L 


8/5/07 


Observed 


163.1986 


57.3838 


1320305-001 


L15 


IRC05 


b;L 


8/5/07 


Observed 






1320306-001 


L15 


IRC05 


b;L 


9/5/07 


Observed 






1320307-001 


L15 


IRC05 


b;L 


9/5/07 


Observed 



in which both changes of filters and dithering in target position 
are performed. The remaining set are for deep surveys and are 
referred to as IRC00 or IRC05, in which no filter change and 
no dithering is performed, in order to minimize the dead time 
and maximize the observation time. For the IRC00 and IRC05 
AOTs, at least three pointed observations are required in order 
to ensure redundant and reliable observations. For these AOTs 
dithering is performed by the operations team by a slight offset 
of each pointing from each other. 

Note that in addition, there exists a jitter between frames in 
the IRC images causing frames to become misaligned with each 
other (note this is not an intentional dither which is a separate 
procedure). The attitude of these frames must be matched (for 
any shift and rotation in position) before stacking to produce the 
final image. 

3.2. Image correction and stacking 

Unfortunately, the FU-HYU data in GOODS-N was taken with- 
out any spacecraft dithering. This means that bad pixels and 



other image artifacts such as slight differences between columns 
in the array, etc always fall in the same positions and are there- 
fore not removed from the final images. In order to overcome 
this we have made use of the intrinsic spacecraft jitter naturally 
present during any AKARI observation to simulate a dithering 
operation. Furthermore, the image distortion stage of the stan- 
dard IRC pipeline toolkit involves many interpolations, with un- 
satisfactory results particularly near bad pixels (see Figure 01. 
Most perniciously, the pixels are no longer statistically indepen- 
dent after this stage. 

To overcome the combined problems of the dither opera- 
tion and the pixel distortion, we customised the standard IRC 
pipeline toolkit to output FITS data products before the field dis- 
tortion stage. At this stage, the data should be dark-subtracted 
and flat-fielded, however bad columns were still evident in some 
of the data in both the L18W and Sll channels. Our approach 
to removing these features exploits the fact that features con- 
fined to single columns cannot be real astronomical signals. 
Initially, after normalising each frame to its median value, we 
fit Gaussian profiles to the histograms of the data values for each 
pixel, and used the inferred modal values to construct a new sky 
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Fig. 3. 100x100 pixel segment of an L18W frame before left and after right the field distortion correction in the standard IRC 
pipeline toolkit. Note the bad pixels which have been spread horizontally over multiple pixels in the right-hand frame, and the subtle 
vertical strobing. 



flat. This technique has been applied in the past to 750-ISOCAM 
data in the ELAIS survey ( jSerjeant eTa l. (2000)) but unlike the 
ELAIS case there is little spacecraft movement between frames. 
Therefore, the flat field constructed in this way will inevitably 
contain real astronomical features as well as the bad columns 
we are seeking to correct. In order to isolate the bad columns 
from the real astronomical point sources we opted to median 
filter the sky flat using periodic boundary conditions, with a 20- 
pixel boxcar window, and subtracted this smoothed image from 
our sky flat. Outlier pixels in the sky flat were replaced with 
global median values for the purpose of this smoothing. The re- 
sulting image represents the fractional excess or deficit of each 
pixel over its neighbours, minus one. We then fit each column of 
this image with a sixth-order Fourier expansion; 



y(x) = Jjiis 



i=0 



i\2n h (/), 

Xmn v 



(1) 



where the row number x runs from 1 to ;t mav =256, the order 
i runs from to 6, and 0, and a, are the i th order phase and am- 
plitude respectively. There are 13 free parameters in total, and 
256 data points in each column. Data points outside the range [- 
0.006,0.006] were excluded which was sufficiently of low-order 
as to not remove point sources from the image, while having a 
sufficiently high order to track the profiles of the bad columns. 
The periodic boundary conditions implicit in the Fourier ex- 
pansion ensured no over-fitting at the image edges. These fits 
gave us a model for the fractional excess flux in each column of 
the detector which was then used to correct the bad columns in 
the FITS data products output from the standard IRC pipeline 
toolkit. 

The next stage of the processing is the determination of the 
intrinsic jitter in the IRC pointed observations. We identified 
known bright sources in the image field, and calculated their 
centroids in every frame. We then corrected the pixel positions 
for field distortion and calculated the mean offset of each frame 
from the first frame, which was assumed to have no jitter offset 
as a reference frame. 

Noise images were then created by estimating the noise 
level in each frame by fitting a Gaussian to the histogram 
of the pixel values. This estimator has previously been 
used in data reduction for ZSO-ISOCAM and at sub-millitre 



Serjeant et al. (2003)1. Pixels greater thanlOcr above the modal 
value or less than -5cr below it were masked, as were their neigh- 
bours and next neighbours. Any regions removed by the standard 
IRC pipeline toolkit were also masked. Masking was achieved 
by assigning an arbitrarily high noise value to the affected pix- 
els. 

As eluded to earlier, the field distortion correction car- 
ried out by the standard IRC pipeline toolkit involves many 
interpolations with every stage of interpolation adding cor- 
related noise to the final products. To avoid such interpola- 
tion steps the IRC pixel fluxes were inserted and coadded di- 
rectly into a zerofootprint map (drawing an analogy to driz- 
zling, c.f. Fru chter & Ho ok(2002)), accounting for the field dis- 
tortion, a procedure again adopted from processing submillime- 



tre data from the SCUBA instrument ((e.g.|Serjeant et al. (2003) 



wavelengths for the SCUBA instrument (Serjeant et al. (2000) 



|Mortier et al. (2005)||Coppin et al.(2006)| l - pictorially shown in 
Figure |U). In total 90 individual L18W frames (one frame was 
discarded due to a moving object in the frame), 180 SI 1 frames 
and 3 1 N4 frames were co-added in this way. 

It was necessary to perform a further stage of deglitching 
to remove negative spikes and we rejected features below -3.5<t 
in the L18W band image and -4cr in the Sll band image. The 
resulting image still showed some subtle large-scale structure, 
possibly a result of poor flat fielding, so we constructed a back- 
ground model by boxcar median smoothing the image with a 
width equivalent to 25 IRC pixels, with periodic boundary con- 
ditions, and subtracted this background model from our images. 
The zerofootprint map was then resampled back to the IRC pixel 
scale. During this process it became evident that the IRC Sll 
data still exhibited structure on smaller scales as shown in Figure 
[5] A possible reason for this is that the IRC pixel scale has 
only 1.68 pixels spanning the expected full-width half maximum 
(FWHM) of a diffraction-limited 1 l//mPSF (Ofwhm = 1-22/1/ D, 
A-l l//m and Z)=0.69m), which is just below the Nyquist scale. 
The same excess structure was not evident in the L18W image, 
therefore we decided to use the rebinned to the original IRC pixel 
scale image at 18/rni for later source extraction but to retain the 
unbinned zerofootprint map at 1 lpm. 

The final co-added images in the AKARI L18W, SI 1 & N4 
bands are shown in the bottom panels of Figure [6] For compar- 
ison in the top panels of Figure [6] we also show the final equiv- 
alent images produced by the standard IRC pipeline processing 
(i.e. for the original undithered observations). The hot pixels that 
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Fig. 4. The methodology for resampling the IRC pixels to a finer 
grid. In the upper panel, the flux of an IRC pixel is inserted into 
a single pixel of a finer grid. Altogether, the 256x256 IRC pixels 
only sparsely populate the finer grid. The lower panel shows a 
later, jittered observation. Due to the jitter, the IRC pixel fluxes 
are not inserted into the same fine pixels as before. As more 
jittered observations are added, the fine grid fills. Where a fine 
grid pixel already contains a measurement, a noise-weighted co- 
addition is made. 



remain after the standard processing are immediately apparent 
in the L18W and Sll band images. Note also the bad columns 
in the standard processed images that have been almost entirely 
corrected for in the new co-added images. 

3.3. Astrometry and source extraction 

Nominally the standard IRC pipeline toolkit calculates astrom- 
etry by matching bright stars in the IRC image to reference ob- 
jects using the 2MASS star catalogue. In reality, this method is 
successful for the IRC-NIR and MIR-S channels, however for 
observations taken with the IRC MR-L channel, the number of 
bright stars is usually insufficient to match to the 2MASS data. 
In this scenario, usually data taken simultaneously in the MIR- 
S/NIR channels is used to project the correct astrometry onto 
the MIR-L images. However, in our case we have already in- 
terrupted the standard IRC pipeline toolkit before the astrome- 
try stage, however FU-HYU by definition has the advantage of 
ancillary multi-wavelength data in place a priori. Therefore, in- 
stead we register the images onto the Spitzer reference frame 
using the brighter 24pm Spitzer sources. 

Source extraction is made using a noise-weighted matched 



filter following the method of Serjeant et al. (2003) For any sig- 
nal image / with associated, noise image N, with point spread 
function P, and a weight image given by W = 1 IN 2 , the mini- 
mum x 2 estimator for a point source flux at any position in the 
images is given by; 



S = 



((/xW)(g)P) 





Fig. 5. Top: A section of the Sll zerofootprint image before 
(left panel) and after (right panel) rebinning to the IRC pixel 
scale. Note the features in the (left panel) at about 9 o'clock 
and 1 o'clock around the two brightest point sources in the im- 
age These features are not apparent in the rebinned image in the 
right panel. Bottom: Stacked point spread function in the Sll 
band, estimated from 101 bright point sources as discussed in the 
text. Contours are spaced at logarithmic intervals of V2 from a 
starting value of 1 / V2. The point spread function is normalised 
to a peak flux of unity. The axes indicate the sizes of IRC pix- 
els, though the image sampling is 4x4 finer. Note the features at 
9 o'clock and 1 o'clock, also visible in the Top left panel , and 
another fainter feature at about 5 o'clock. Companion galaxies 
or noise spikes may cause the weak features in the bottom-right 
corner and top-centre of the image. 

where (^) denotes a convolution. The associated error on this 
quantity is given by; 



AS = 



1 



(3) 



,p2 



,p2 



(2) 



The ll//m PSF was estimated directly from the image, us- 
ing 101 bright point sources in the zerofootprint image (see bot- 
tom panel of figure [5]). Inevitably, some of the point sources 
used to compile the PSF will have neighbouring objects. To 
alleviate this problem, we used an algorithm originally devel- 
oped for estimating mean quasar fluxes with nearby compan- 
ions (Serjeant & Hatziminaoglou (2009)). We assume our mea- 
surements Xi with estimated noise levels cr, have an additional 
unknown noise contribution cr„ caused in this case by the RMS 
fluctuations from companion galaxies. The maximum-likelihood 
solution is then determined for the mean level p and for cr si- 
multaneously. When using this technique on all observations of 
any given pixel of the point spread function, it is found to be a 
less noisy estimator than for example, the median value. We re- 
peated this technique for all pixels of the point spread function. 
While this procedure is less sensitive to outliers caused by single 
companions, it remains possible that some excess flux remains, 
nevertheless, if such a low-level excess is typical around our ob- 
jects, it is appropriate and indeed necessary to include this in the 
point source filter. We neglect any possible variation in the PSF 
over the field of view at the present time. Since the L18W im- 
age does not suffer from the same small scale structure we can 
assume a Gaussian PSF with the standard instrument estimation 
of the FWHM equal to 2.3 pixels. 
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Fig. 6. Top: The final processed images produced by the standard IRC pipeline; (left panel) L 1 8 W band image, (middle panel) S 1 1 
band image, (right panel) N4 band image. Bottom: The final processed images produced by terminating the IRC pipeline before 
the distortion correction step and using the spacecraft jitter to provide dithering to correct for hot pixels described in the text; (left 
panel) L18W band image, (middle panel) SI 1 band image, (right panel) N4 band image. Note the absence of hot pixels (appearing 
as black spots in the Top panels in the standard processed images) particularly in the L18W and SI 1 band images. The banding and 
bad columns in the original images are also greatly improved. 



The best fit flux at any point in the image map is given by 
the value in the PSF convolved image at the extracted source po- 
sition and the corresponding errors from the values in the noise 
map. Sources are extracted by thresholding the signal-to-noise 
(o~ - S/AS) image. In practice, objects that are separated at 
a high signal-to-noise threshold may be connected at a lower 
threshold. We therefore thresholded the signal-to-noise images 
from 3cr to 10cr in cr steps, using the locations of the peak signal- 
to-noise as the source position, and selected only the unique 
sources from the resulting lists to produce our source catalogues. 
The final signal-to-noise images, before and after PSF convolu- 
tion, are shown in Figure [7] 

Note that although the standard AKARI IRC Instrument man- 
ual ( |Lorente et al. (2007)) > provides conversion factors from the 
raw instrument ADU units to Jy, this conversion factor is based 
on aperture photometry of standard stars with a flux aperture and 
assumed sky annulus. This conversion factor is not appropriate 
for our source extraction algorithm and we have thus returned to 
the original calibration stars (see Tanab e et al. (20 08) ) and repro- 
cessed the raw calibration data using exactly the same method as 
applied to the GOODS-N data in order to calculate the appropri- 
ate conversion factors for our photometry method. 

In our final extracted source lists we find 233 sources at >3cr 
in the L18W image with fluxes >100^/Jy, and 242 >3cr sources 
in the Sll image with fluxes >50^Jy. The L18W survey covers 
around 101.3 arcmin 2 with l<x noise levels below 100/Jy, and 
98.3 arcmin 2 below 50/iJy. The Sll ll/mi survey covers 98.06 
arcmin 2 with noise levels below 100/Jy, and 95.6 arcmin 2 below 
50^/Jy. The median noise levels at 18/vm and ll/vm are 37.7/Jy 



and 23.0yuJy respectively. These figures assume conversions of 
1 ADU/s = 27.36/Jy at 18/mi and 22.23/Jy at 11p.m. There is 
a possible absence of sources in the ll/zm and 18pm images, 
tracing the bad columns. It is possible that an over-zealous bad 
column correction has removed astronomical flux. However this 
is also seen in reduction using the standard IRC pipeline toolkit. 
The causes (assuming it is an artefact) are still under investiga- 
tion. In the N4 band we detect 340 sources down to a 3cr flux 
limit of ~8yuJy. 

4. Cross associations and catalogue creation 

4.1. Cross-associations with Spitzer IRAC & MIPS data 



As part of the FU-HYU Mission program, Negrello et al. (2009) 
have cross-associated the Spitzer data for the GOODS-N field 
with our AKARI catalogue. The Spitzer Legacy Data Product 
release 3 are available in the public domairQ and consist of 
both the images in the four IRAC bands (3.6, 4.5, 5.8, 8.0pm) 
and the image and source catalogue from the MIPS 24pm 
band. Sources were extracted from the four IRAC images us- 
ing SExtractor (Ber tin & Arnouts(1996)| l resulting in catalogues 
containing 5,792, 5,576, 2,328, 2,186 3cr sources in the 3.6, 
4.5, 5.8 & 8.0/im bands respectively. Sources with SExtractor 
'STAR.CLASS' equal to eit her 1 (cosmicray-like) or (star- 
like) were then removed (See Negrello et al. (2009) for details). 
The MIPS source catalogue consists of 1,199 3<r sources to a 
flux limit of 80 piy. However, the MIPS24 image obviously in- 



http://data.spitzer.caltech.edu/popular/goods/Documents/goodsdataproducts 
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Fig. 7. Final GOODS-N signal-to-noise images in the L18W band (left panels), Sll band (middle panels) & in the N4 band (right 
panels) , before the point source filtering (top, scaled from -5<x to 10<x) and after (bottom, scaled from -10cr to 20cr). 



eludes more sources than the catalogue and we have reprocessed 
the image using DAOPHOT to extract sources deeper than the 
80 yuJy limit. 

As a reference for our band merged catalogue we have 
started with the catalogue created from the deep IRAC 3.6/im 
band image and cross-correlated all catalogues with this. 
Sources between catalogues were matched over a radius of 

r = ^(o~\ blim + cr 2 A ), where <x is twice the Gaussian rms width 
of the instrument beam (equivalent to the Full Width at Half 
Maximum, FWHM/(2 V2Zn2)) and A is the wavelength of the 
cross-match catalogue. Only 3.6/im sources with a single unique 
counterpart in at least one of the other catalogues were included 
in the band-merged catalogue. When no counterparts were found 
in a given catalogue the flux and error were set to -1. Finally any 
pairs of objects sharing the same flux were removed from the 
final band-merged catalogue. 

The final band merged catalogue contains a total of 4393 
sources at 3.6pm. There are 4344, 2045 & 1873 associations in 
the 4.5, 5.8 & 8pm IRAC bands respectively. In the MIPS 24pm 
band we find 1 143 associations. In the AKARI bands we find 
340, 181 & 192 associations in the N4, Sll & L18W bands re- 
spectively. 

For the AKARI sources in the Sll band 13 llyum sources 
and in the L18W band, 12 18^m sources have no counterpart 
at 3.6^m respectively and are therefore not included in the cat- 
alogue. On visual inspection we find that almost all of these 
sources either lie on bright sources or around the periphery of 
the AKARI image making their reliability dubious. 

In Table [3] the sources with detections in all 3 AKARI bands 
are listed with fluxes in piy for the Spitzer IRAC 3.6, 4.5, 5.8, 
8jum, IRS 16/an, MIPS24 bands, the AKARI IRC N4, SI 1, L18W 
bands and the ISO 15pm band. The full FU-HYU catalogue will 
be made available to the astronomical community in the future. 



4.2. Cross-associations with other data 

In addition to the Spitzer IRAC & MIPS data, the GOODS- 
N field enjoys copious multi-wavelength coverage at other 
wavelengths. Using the same matching criteria described above 
we have added additional bands to our FU-HYU-GOODS- 
N catalogue. A 36 square arcminute area of the GOODS-N 
field was also imaged by Spitzer at 16pm down to a flux 
level of ~0.08mJy by Teplitz et al. (2005) using the Infra-Red 



Spectrograph (IRS) blue peak-up filter, detecting around 153 
sources. Cross-correlating with our data set we find a total of 83 
sources in our catalogue. The GOODS-N field was also imaged 
by ISO at 15pm, down to 0.2mJy detecting of the order of 100 
objects ( Ausse let al.(1999")) . We find numerous close pairs with 
similar flux levels and after selecting unique sources we find 59 
associations with our catalogue. 

The GOODS-N field has also been imaged at optical 
wavelengths in the U,B,V,R,I,z' bands over the entire field 
and partially in H,K' using both space-bourne Hubble Space 
Telescope / Advanced Camera for Surveys (HST/ACS) and 
the Subaru, KPNO & Hawaii 2.2m ground based observato- 
ries ( |Capak et al.(20 04) ). we find approximately 3600 sources 
in the UBVRIz' bands associated with sources in our catalogue. 
Specifiacally, at the UBVRI, z', HK' detection limits of 27.1, 
26.9, 26.8, 26.6, 25.6, 25.4, 22.1, 22.1 AB magnitudes we asso- 
ciate 3536, 3609, 3617, 3586, 3617, 3612, 3380, 3380 sources 
in our catalogue. 

In the ultra-violet, the GOODS field has been observed as 
part of the surveys made by the Galaxy Evolution Explorer 
(GALEX, |Martin et al. (2005)) in two bands (232nm and 
154nm) referred to as the near-uv (nuv) and far-UV (fuv). For 
cross-associations with our FU-HYU catalogue we use the data 



of Morrissey et al. (2007) with fluxes derived from the data re- 
duction pipeline of Burgarell a et al.(2007)| In total we find 829 
associations with our catalogue in both the GALEX bands. 
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Table 3. Source catalogue in the FU-HYU GOODS-N survey for sources detected in all IRC bands with fluxes given for the Spitzer 
IRAC 3.6, 4.5, 5.8, 8pm, IRS 16pm, MIPS24 bands, the AKARI IRC N4, SI 1, L18W bands and the ISO 15pm band. 



ID 


RA 


DEC 


IRAC3.6 


IRAC4.5 


IRAC5.8 


IRAC 8 


IRS 16 


MIPS24 


N4 


Sll 


L18W 


IS015 




(deg.) 


(deg.) 


0"Jy) 


(pJy) 


(pJy) 


G"Jy) 


G"Jy) 


(pJy) 


0"Jy) 


(pJy) 


0"Jy) 


(pJy) 


357 


189.233 


62.136 


108.50 


76.00 


73.21 


73.86 


568.00 


798.05 


82.17 


213.41 


383.99 




417 


189.194 


62.143 


111.38 


80.98 


72.10 


84.81 




990.90 


104.39 


115.08 


599.85 




447 


189.138 


62.143 


102.71 


90.25 


91.81 


109.02 




769.89 


84.15 


164.61 


428.38 




456 


189.256 


62.145 


37.67 


20.72 


23.05 


21.29 


130.00 


184.85 


24.03 


69.56 


122.91 




481 


189.209 


62.146 


101.31 


86.28 


79.78 


129.94 


348.00 


561.86 


80.18 


323.41 


292.51 




485 


189.151 


62.145 


50.99 


36.53 


36.47 


35.01 




189.86 


36.91 


179.11 


148.63 




500 


189.167 


62.145 


135.39 


86.98 


82.80 


426.77 




313.06 


80.06 


326.84 


266.81 




533 


189.127 


62.148 


77.82 


52.27 


48.90 


40.32 




380.92 


48.77 


99.75 


254.62 




884 


189.320 


62.169 


75.16 


66.17 


54.47 


132.94 




579.73 


81.09 


361.80 


299.22 




1082 


189.121 


62.180 


98.14 


88.33 


76.03 


70.62 


465.00 


720.49 


84.08 


75.10 


433.77 




1101 


189.349 


62.180 


92.92 


62.13 


54.92 


242.26 


141.00 


182.99 


60.26 


183.92 


128.81 




1261 


189.089 


62.186 


136.10 


102.12 


80.44 


68.78 






123.53 


83.09 


290.68 




1355 


189.228 


62.191 


50.10 


43.73 


27.11 


118.23 




166.28 


35.55 


180.29 


103.31 


42.00 


1381 


189.330 


62.192 


74.11 


52.31 


48.52 


50.71 


213.00 


201.20 


54.04 


171.39 


200.44 




1441 


189.078 


62.198 


32.97 


22.06 


22.24 


14.47 




155.75 


27.19 


60.66 


99.02 




1464 


189.154 


62.193 


170.63 


114.91 


113.53 


545.24 




788.55 


94.52 


370.46 


376.82 


300.00 


1471 


189.257 


62.196 


60.04 


70.31 


69.05 


54.47 




673.02 


59.82 


65.06 


384.33 


15.00 


1590 


189.246 


62.203 


62.69 


42.63 


35.08 


24.90 




264.50 


45.69 


85.11 


179.07 


157.00 


1680 


189.316 


62.200 


300.51 


192.04 


126.56 


101.12 




217.99 


210.98 


89.12 


99.20 




1698 


189.090 


62.208 


47.05 


42.06 


31.48 


72.24 




251.90 


47.72 


199.56 


169.25 




1716 


189.153 


62.204 


126.41 


83.66 


68.90 


50.37 


300.00 


395.21 


88.33 


101.26 


260.65 


202.00 


1760 


189.066 


62.210 


70.31 


64.61 


41.47 


155.21 




394.12 


57.64 


263.27 


232.65 




1771 


189.144 


62.204 


316.28 


244.52 


196.02 


338.39 


853.00 


1309.87 


227.87 


858.14 


896.83 


448.00 


1827 


189.225 


62.215 


58.38 


37.31 


36.45 


25.33 


207.00 


134.91 


38.75 


91.85 


180.88 


179.00 


1885 


189.090 


62.217 


41.04 


28.58 


31.71 


27.42 






35.14 


61.65 


129.56 




1907 


189.081 


62.215 


129.55 


105.91 


100.16 


167.43 




975.50 


106.52 


467.30 


595.12 




1985 


189.140 


62.222 


49.92 


32.84 


34.63 


23.08 


320.00 


319.17 


40.65 


68.49 


192.08 


122.00 


2109 


189.096 


62.230 


39.30 


35.63 


30.37 


69.09 




296.78 


37.49 


236.64 


247.21 




2253 


189.074 


62.236 


52.69 


57.54 


56.07 


46.33 




447.95 


58.73 


61.34 


300.72 




2345 


189.041 


62.240 


63.87 


40.00 


42.13 


26.40 




293.24 


48.45 


91.51 


199.19 




2355 


189.058 


62.238 


59.06 


48.61 


36.09 


213.25 




222.45 


55.01 


260.37 


215.32 




2412 


189.201 


62.241 


86.26 


60.13 


54.77 


386.56 


283.00 


463.48 


68.03 


363.90 


290.78 


307.00 


2483 


189.148 


62.240 


71.64 


99.93 


163.45 


282.56 


615.00 


1487.41 


89.34 


368.62 


755.66 


441.00 


2558 


189.250 


62.247 


50.47 


33.56 


39.29 


32.95 




453.35 


37.32 


165.57 


310.12 


295.00 


2769 


189.096 


62.257 


46.39 


51.08 


72.85 


130.06 


335.00 


519.55 


50.95 


159.34 


347.93 




2788 


189.165 


62.257 


46.59 


40.28 


27.90 


66.68 


121.00 


139.45 


39.61 


120.86 


135.86 




2885 


189.193 


62.258 


84.54 


56.75 


54.15 


42.26 


433.00 


524.84 


59.16 


119.42 


350.03 


418.00 


2955 


189.132 


62.268 


44.94 


29.99 


30.94 


23.77 


245.00 


286.83 


32.54 


86.56 


238.93 




2979 


189.261 


62.262 


143.18 


131.91 


132.58 


158.66 




482.65 


125.25 


265.34 


297.76 




3023 


189.176 


62.263 


126.59 


102.57 


106.10 


120.96 




829.52 


129.85 


190.16 


457.31 


459.00 


3073 


189.252 


62.271 


41.08 


29.99 


32.23 


28.41 




490.26 


31.41 


78.12 


285.35 




3107 


189.145 


62.275 


58.86 


40.10 


38.46 


34.78 


368.00 


458.81 


41.26 


104.18 


299.68 




3191 


189.245 


62.277 


81.26 


70.72 


48.76 


148.61 




282.12 


64.82 


272.17 


260.71 




3233 


189.158 


62.297 


162.41 


120.65 


72.01 


132.01 




385.86 


108.53 


294.75 


313.63 




3398 


189.225 


62.294 


80.09 


72.90 


56.07 


143.53 




349.00 


72.39 


345.16 


318.29 




3479 


189.183 


62.285 


30.85 


25.36 


20.70 


42.54 




210.30 


19.17 


95.30 


114.45 




3725 


189.277 


62.292 


94.69 


83.88 


57.85 


150.07 




303.69 


76.49 


285.99 


242.33 





5. Results 

5.1. Final catalogue 

The final catalogue contains data for 4393 sources over 19 bands 
from mid-infrared to UV wavelengths and their associated er- 
rors. However, in the final product, we find only 4 sources with 
data in all 19 bands. Concentrating on our near to mid-infrared 
sources only we find a total of 7 sources with fluxes in all 10 
infrared bands (i.e. Spitzer IRAC 3.6, 4.5, 5.8, 8pm, IRS 16pm, 
MIPS 24pm, AKARI IRC N4, Sll, L18W & ISO 15pm bands). 
Two of the 3/7 sources without data in all 19 bands have no de- 



tection in the GALEX bands and the other of the 3/7 sources has 
no detection in either the GALEX bands or optical wavelengths. 

The total number of sources with complete infrared coverage 
increases to 17 sources if the constraint on the ISO detection 
is relaxed. Removing both the relatively smaller 16pm Spitzer- 
IRS and 15pm ISO samples we retrieve a sample of around 50 
sources. Whilst relaxing the constraint on the AKARI N4 band 
we are left with a sample of 100 sources with complete coverage 
of their infrared spectrum from 3-24pm (See Figure[TJi. We refer 
to this hereafter as the the total mid-infrared (MIR) flux (the total 
flux in the Spitzer IRAC 3.6, 4.5, 5.8, 8pm, AKARI Sll, L18W 
& Spitzer MIPS 24pm bands). 
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Fig. 8. Infrared colour-colour diagram for sources n the FU- 
HYU-GOODS-N field with and without optical counterparts. 
Colours are IRAC 3.6^m / IRC L18W & IRC L18W / IRAC 
5.8/im. There is no discernible difference in the colour-colour 
space populated by the sources without optical counterparts. 



There is also a sample of 27 sources with no optical coun- 
terpart in any of the optical bands with Spitzer IRAC, MIPS and 
AKARI MIR-L identifications. To investigate the significance of 
this population we plot in Figure [8] the infrared colours of the 
sources without optical counterparts along with the sources with 
optical counterparts. It is found that the two populations occupy 
the same colour-colour space with no discernible differences. 

5.2. Sources with full infrared coverage 

The spectral energy distributions (SED) of the 7 sources with 
complete infrared coverage in all 10 AKARI, Spitzer & ISO in- 
frared bands are shown in Figure [9] and listed below. In Figure 
[9] the best spectral fit and photometric redshift from a com- 
bination of starburst and AGN templates using the algorithm 
of Negrello et al. (2009) is shown. The starburst templates are 
taken from the libraries of Takagi et al.(2003)| which use three 
types of dust model to model the extinction, referred to as 
the Milky Way (MW), Large Magellanic Cloud (LMC), and 



Small Magellanic Cloud (SMC) dust models and are a func- 
tion of the metallicity. Model parameters are the star forma- 
tion timescale (lOOMyr), metallicity (assumed to be 0.1), the 
compactness of the starburst region (0) and the starburst age. 
The variation in the spectra is explained by the difference in 
the age (ranging from 0.01-6Gyr) and the compactness of the 
starburst region (ranging from 0.3 - 5.0 which can be consid- 
ered as a measure of the optical depth with values of of 
0.3 and 5.0 corresponding to high and low optical depths re- 
spectively). The AGN templates are taken from the libraries 
of Efstath iou & Rowan-Robinson (1995)| who assume a tapered 
disc geometry for the AGN with the the torus height increas- 
ing with distance from the central source, tapering to a constant 
value at large distances. The tapered disc is assumed to have a 
1 jr density distribution and a ratio between the outer and the in- 
ner radius of the torus of 20. The spectra of the torus depends 
only upon the viewing angle of the torus, v i ew , with # v iew=0°, 
90° corresponding to an edge on and face on torus respectively. 
The model parameters used for the fits, and the resulting redshift 
and luminosity contributed by the starburst and AGN compo- 



nents respectively are tabulated in Table [4] All the quoted errors 
represent the 99 per cent confidence interval. 

FU-HYU ID=G1716 This source has a 24/zm flux of 395±8yuJy 
and an ISO 1 5/rni flux of 202+58 /Jy (HDF PM3 8 in the 
source list of Aussel et al.(1999)). The best spectral fit is 
a PAH dominated starburst (i.e ULIRG with L~10 12 L©) at 
redshift 1.02. 

FU-HYU ID=G1771 This is a bright 24pm source with flux of 
1310+12juJy. The IS O 15/rni flux is 448+ 68/Jy (HDF PM3 
2 in the source list of |Aussel et al.(1999)| ). The best spectral 
fit is a PAH starburst galaxy (i.e ULIRG with L=10 12 3 L o ) 
with underlying mid-infrared AGN at redshift 0.78. 

FU-HYU ID=G1827 This source has a relatively faint 24pm 
flux of 134.91+6yuJy. The ISO 1 5/zm flux is 179+6 0/Jy 
(HDF PM3 33 in the source list of |Aussel et al.(1999)| >. The 
best spectral fit is a PAH dominated luminous infrared galaxy 
(LIRG, L=10 n A L Q ) at redshift 0.76. 

FU-HYU ID=G1985 This source has a 24^111 flux of 
319.7il2yU.Jy. The ISO \5 pm flux is 122+54 /Jy (HDF 
PS 3 3 in the source list of | Aussel et al.(1999)) . The best 
spectral fit is a pure PAH dominated luminous infrared 
galaxy (LIRG, L=10 115 L o ) at redshift 0.88 with strong 
silicate absorption. 

FU-HYU ID=G2412 This source has a 24^111 flux of 
463+7juJy. The ISO \5p m flux is 307+62/ Jy (HDF 
PM3 24 in the source list of |Aussel et al.(1999)| >. The best 
spectral fit is a pure PAH dominated moderate starburst 
(L=10 101 L o ) at redshift 0.12. 

FU-HYU ID=G2483 This source is the brightest infrared 
source in our sample of complete coverage with a power 
law spectrum and 24pm flux of 1487+1 1//Jy. The ISO 
\5pm fl ux is 441+43/Jy (HDF PM3 5 in the source 
list of |Ausseletal.(1999)| ).The best spectral fit from 



Negrello et al. (2009) is an AGN dominated source with a 



featureless mid-infrared pectrum at redshift 0.04. Although 
the model fitting predicts and roughly equal partition be- 
tween the starburst and AGN luminosities, the uncertainty 
on the starburst luminosity is large. This source has no de- 
tections in the GALEX UV bands. 
FU-HYU ID=G2885 This source has a 24pm flux of 
524±8//Jy. The ISO \5pm fl ux is 418+9huJy (HDF 
PM3 21 in the sou rce list of | Aussel et al.( 19 99) ). The 



best spectral fit from Negrello et al. (2009) is a composite 



starburst & AGN at redshift 1 .0 and would be classed as a 
ULIRG with L=10 12 Lq . This source has no detections in 
either the UBVRIz'HK' optical bands or the GALEX UV 
bands. 



5.3. Source fluxes and colours as a function of the total 
mid-infrared emission 

Given the comprehensive multi-wavelength mid-infrared cover- 
age provided by the combined AKARI & Spitzer bands, we de- 
fine the quantity referred to as the total mid-infrared flux as the 
total flux in the Spitzer IRAC 3.6, 4.5, 5.8, 8pm, AKARI SI 1, 
L18W & Spitzer MIPS 24pm bands for the sources in the FU- 
HYU-GOODS-N field. This quantity is representative of the to- 
tal emission from the mid-infrared spectra, dominated by the 
PAH emission features in our sources. In Figure [TO] this total 
mid-infrared flux is plotted as a function of the mid-infrared 
fluxes (panel (a)) and colours (panel (b)) in the Spitzer-IRAC 
8pm & -MIPS 24pm bands and the AKARI-IRC Sll, L18W 
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Table 4. Best fit results for the 7 objetcs in Figure[9] 
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log[L„ e „] 
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(deg) 
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(AD 
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n n+0.16 



1771 0.78^f 4.45 0.00 13 SMC 400 + || l.6 + }\ 42 + l 12.3+^ 10.9 + | : | 12.3 + |^ 

1827 0.76 + U :i )| 2.65 0.00 13 LMC 500+]™ 1.6 + ° : j 47+f 11. 4 + J™ 9.92 + ||* H-4 + ™ 

1985 0.88+° : ^ 1.33 0.19 11 MW 500.!™- 2.0™ 50« 2 11.5+°g 9.87^™ 1L5 -0Q6 

2412 0.12+° : * 1.35 0.17 13 MW 600!%' 3.0^ : | 47+| 10.(T° : g 8.52+^ 10.1^' 

2483 0.04^ 1.40 0.16 11 LMC 10^ 90 1J2«| 45« 9.67^ : ^ 10.7! 2: " 

2885 0.87 0.47 4 SMC 600!™ 0.7+^ 8!f LL5^| 12-0!g;jg 

Zphot^ Best fit photometric redshift 

/tmfa/ v: reduced minimum^ 2 for v degrees of freedom 

f„2: probability associated to the minimum^ 2 

ext.: extinction curve 

Age: age of the starburst 

0: starburst compactness factor 

$ v j cw : viewing angle of the AGN torus 

log[L s j,]: Luminosity contributed by the starburst component 
\og[L agn \: Luminosity contributed by the AGN component 
log[L, D ,]: Total Luminosity 
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Fig. 9. Photometric data for the seven sources in the FU-HYU GOODS-N field with fluxes in all 10 infrared bands (i.e. Spitzer 
IPvAC 3.6, 4.5, 5.8, 8/im, IRS 16//m, MIPS 24pm, AKARI IRC N4, S 1 1 , LI 8W & ISO 15/rni bands). Also shown are the model fits 
for these sources using the photo-metric redshift code of Negrello et aL (200 9)| The dashed and dotted lines are th e fitted starburst 
and AGN components from the spectral libraries of Takagi et al.(2003) & Efstathiou & Rowan-Robinson (1995) respectively and 
the thick solid line is the composite fit. Photometric redshift fitting parameters are listed in each panel and tabulated in Table |4] 
following [Negrello et al. (2009)| 



bands respectively. We observe a tight correlation between the 
total mid-infrared flux and the longer wavelength MIPS 24/mi 
and IRC L18W bands increasing as a function of increasing to- 
tal MIR flux but a larger scatter in the shorter wavelength bands 
revealing that the longer mid-infrared bands are better indica- 
tors of the total mid-infrared flux of these sources. The colours 
exhibit an almost constant value for the the case of L18W/MIR 
& MIPS24/MIR as a function of total MIR flux where as the 
IRAC8/MIR, S 1 1/MIR colours show a slight redder trend to high 
values of the total mid-infrared flux. The corresponding near- 
infrared flux is plotted in Figure QT| as a function of the mid- 
infrared fluxes (panel (a)) and colours (panel (b)) for the Spitzer- 
IRAC 3.6,4.5, 5.8pm & AKARI-IRC N4 bands. In general there 
is less dispersion in both the fluxes and the colours (for all IRAC 



and the IRC-N4 bands) as a function of total mid-infrared flux 
since the near-infrared spectra of these sources are less affected 
by the emission features. In general there is a trend of increas- 
ing near-infrared flux as a function of total mid-infrared flux but 
bluer colours in the colour-total MIR flux plane. 

In Figure [T2]the total mid-infrared flux is plotted against the 
optical fluxes (panel (a)) and colours (panel (b)) for the sources 
in our catalogue. For clarity we include only 3 bands (U, R, z')- 
A wide dispersion in optical fluxes is seen over the entire range 
of total MIR flux implying that the two are not particularly cor- 
related in general. There is a trend to increasing optical flux as a 
function of increasing total MIR emission although as expected, 
it is not as pronounced as in the infrared band case. Similarly, 
there appears little or no correlation of optical/MIR colour as a 
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Fig. 16. Model spectral templates adopted for the colour tracks 
plotted in Figures [T7]& [18] and described in the text. 



function of total MIR flux, in particular the U-band/MIR colour 
is distributed over 2 orders of magnitude at all values of the total 
MIR flux. Figure [13] shows the corresponding flux (panel (a)) 
and colour (panel (b)) for the GALEX UV bands. Similarly to 
the U-band plots, no significant correlation is seen in either the 
flux distribution or colours as a function of total MIR flux, al- 
though the near-UV (NUV) emission shows less dispersion than 
the far-UV (FUV) emission. 



5.4. The nature of the galaxy populations in the 
FU-HYU-GOODS sample 

In Figures [14] & [15] we plot a selection of colour-colour distri- 
butions from the FU-HYU-GOODS sample covering the wave- 
length range from 3.6-24yum in the Spitzer-IRAC 3.6, 8pm 
bands, the AKARI-IRC SI 1, L18W bands and the Spitzer-MlPS 
24pm band. Figures \14\a) & (b) plot the above fluxes as a func- 
tion of^wjc/IRAC8 - Total MIR/ flux colour andfluxflRC Sll 
- Total MIR/ flux colour respectively. A distinct separation is 
seen between the longer L18W & MIPS24 band colours (which 
show very similar colours and trends) and the shorter wave- 
length IRAC3.6 & IRC SI 1 band colours. The anomalous L18W 
& MIPS24 band point at lg(MIR/L18W) & lg(MIR/MIPS24) 
colour ~ 1 .7 is a bright local galaxy. Note that although the IRAC 
band colours occupy similar colour-colour space, the AKARI- 
IRC SI 1/IRAC8 colours in Figure \\A\g) extend to redder values 
(bluer values on the Total MIR/S 1 1 axis) implying an excess in 
the Sll band emission not seen in the shorter wavelength IRAC 
bands. 

In Figures [TBTa) & (b) the above fluxes are plotted as a 
function of the longer wavelength^Mx/MIPS24 - Total MIR/ flux 
colour andfluxflRC LI 8W - Total MIR/ flux colour respectively. 
Similarly to Figure [14] the same separation in colour-colour 
space is seen for the longer (MIPS 24, IRC L18W) wavebands 
from the shorter wavelength IRAC bands. The tight clustering 
of the MIPS24/L18W - Total MIR/MIPS24 in FigureESM and 
the L18W/MIPS24 - Total MIR/L18W in Figure \T5\b) empha- 
sis the fact that the longer wavelength bands provide an accurate 
measure of the Total MIR emission. In contrast the shorter wave- 
length bands show a wider dispersion as a function of colour 
and indeed in some cases exhibit a bi-modal distribution. Note, 



the apparent excess in the Sll band flux seen in Figure \\A\a) 
manifests itself by producing two distinct populations in the 
S 1 1/MIPS24 - Total MIR/S 1 1 and S 1 1/L 1 8 W - Total MIR/S 1 1 
planes, one population with similar colours to the IRAC band 
colours and another population with colours more akin to the 
longer wavelength (MIPS24, L18W) band colours. 

In Figure [17] a selection of colour-colour distributions are 
individually extracted from Figures [14] & Q3] Overlaid on each 
individual colour-colour distribution are the galaxy template 
spectral energy distributions for a normal quiescent galaxy, a 
star-forming galaxy, an ultraluminous infrared galaxy and an 
AGN. The normal spectral template is taken from the SED li- 
braries of Efstathio u & Rowan-Robinson (2003)] which as sume 
an input radiation field based on Bruzual & Chariot (1993)] and 
assumed interstellar dust field. There are two free parame- 
ters, the ratio of the radiation field to the local solar neigh- 
borhood, X, an d the visual extinction Ay. Our adopted tem- 
plate corresponds to ^=5 , Ay=0.5. The star-forming galaxy 
and ultraluminous infrared galaxy templates are taken from 
the libraries of Efstathiou et al. (2000) The templates assume 
a Salpeter initial mass function (IMF) from 0.1-125Mq with 
an e-folding time of 20Myr. The evolution of the stellar pop- 
ulation within giant molecular clouds follows the stellar syn- 
thesis population models of Bruz ual & Chariot (1993)| The in- 
dividual starburst models are defined by two free parameters: 
the age of the starburst in Myr; Tage and the initial opti- 
cal depth (tv) of the molecular clouds. We select two spe- 
cific templates that represent the archetypical starforming galaxy 
M82 and ULIRG Arp 220 which correspond to model pa- 
rameters of Tage=10 & 26 Myr and Ty-50 & 150 respec- 
tively. For M82 we replace the model SED with the 18-5/vm 
region measured by the ISOCAM Circular Variable Filter and 
the true optical SED ( |Forster Schreiber et al. (2003)). Finally 
our AGN SED is derived from the dust torus models of 
Efstathio u & Rowan-Robinson (1995)| with the near- infrared 



optical spectrum following King & Rowan-Robinson (2003) 



The dust torus model assumes an opening angle of 45°, ratio 
of outer to inner torus radius of 20 and a viewing angle of 10° . 
The adopted model spectral templates are shown in Figure[T6l. 

The SED tracks in Figure [T7] are marked by a large 
square at the redshift =0 position and have crosses for ev- 
ery 0.2 steps in redshift thereafter extending to a redshift of 
2. Figure Wf\a) shows the IRAC 8yum/L18W against the Total 
MIR flux/ IRAC 8/im colours. The shaded area bounded by 
lg(IRAC8/L18W)<-0.5, lg(MIR/IRAC8)>1.2 denotes sources 
that are expected to be higher redshift (z>0.7) star-forming 
galaxies. These colours can be compared with similar work 
carried out by |Brand et al. (2 006) using IRAC8, MIPS24 band 
Spitzer data. Our colours agree well with the corresponding seg- 
regation of low and high redshift star-forming population to red- 
der and bluer colours respectively either side of the locus of 
AGN colours lg(IRAC8/L18W)<-0.3. In Figure \ITib) the IRC 
Sll / L18W colour is plotted against the total MIR / IRC Sll 
band colour. A similar trend can be seen in these colours as 
Figure W7\a ) and the colours produce an effective means of seg- 
regating a star-forming population from the normal galaxy pop- 
ulation at high (lg(MIR flux / Sll flux)>0.1, z>l) redshifts. 
Note that this result can explain the bimodal distribution seen 
in FigurdT5lfor the S11/L18W, Sll/MIPS colours, as a popula- 
tion of lower and higher redshift source populations. The multi- 
band coverage of the mid-infrared spectra of our sources can 
also be used to provide powerful diagnostics of the PAH emis- 
sion in our galaxy spectra using colour-colour information. In 
Figure flTTc) the IRAC 8//m / IRC Sll band colour is plotted 
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(b) MIR Colours as function of 
total mid-infrared emission 
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Fig. 10. (a): Mid-infrared source fluxes in the Spitzer & AKARI bands as a function of the total mid-infrared flux (the total flux in 
the Spitzer IRAC 3.6, 4.5, 5.8, 8pm, AKARI Sll, L18W & Spitzer MIPS 24pm bands) for the sources in the FU-HYU-GOODS- 
N field, (b): Mid-infrared colours (mid-IR flux/Total MIR flux) as a function of the total mid-infrared flux for the sources in the 
FU-HYU-GOODS-N field. 
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(b) NIR flux as function of 
total mid-infrared emission 
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(b) NIR colours as function of 
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Fig. 11. (a): Near-infrared source fluxes in the Spitzer & AKARI bands as a function of the total mid-infrared flux (the total flux in 
the Spitzer IRAC 3.6, 4.5, 5.8, 8yum, AKARI Sll, L18W & Spitzer MIPS 24/zm bands) for the sources in the FU-HYU-GOODS-N 
field. fluxes, (b): Near-infrared colours (near-IR flux/Total MIR flux) as a function of the total mid-infrared flux for the sources in 
the FU-HYU-GOODS-N field. 



against the total MIR / IRAC 8/mi colour. The two highlighted 
regions on the plot show the passage of prominent (predomi- 
nantly 1.1pm) PAH features through the AKARI Sll band. From 
redshifts of ~0.4 -0.6 the 1.1pm feature enters the Sll band 
and a sharp decrease in the IRAC 8pm / IRC Sll colour is ob- 
served in the SEDs of starforming (starburst + ULIRG) galax- 
ies. At redshifts > 0.8-1 the feature is redshifted out of the 
Sll band and a gradual climb back up the IRAC 8pm / IRC 
Sll colour track is seen in these sources. Note that for clar- 
ity, the AGN track is not shown on this plot although the tracks 
are constrained to lg(IRAC8/Sll)>-0.15 & lg(MIR/IRAC8)<l. 
Finally, in Figure \VRd) the MIPS24/L18W colour is plotted 
against the total MIR/MIPS24 band colour. Interestingly, it is 
found that following our SED templates, the quiescent normal 



galaxies can be well separated from the evolving galaxy pop- 
ulation by their MIR/MIPS24 colour alone. This could be in- 
dicative of the 24pm flux as an excellent tracer of star forma- 
tion in galaxies where the 24pm emission in normal quiescent 
sources would be systematically lower with a greater dispersion. 
Note that Bav ouzet et al. (200 8) have indicated that the 24pm 
flux is a good indicator of the mid-infrared and following this the 
bolometric infrared luminosity of star forming galaxies (STFG) 
where as |Dale et al.(2007)| have shown for a limited sample of 
local quiescent galaxies a larger dispersion in the mid- to far- 
infrared flux ratios. In addition, for the star-forming populations 
in Figure [TTl't/J a blue tail to the MIPS24/L18W colours is ob- 
served. This is sparsely populated and almost certainly due to 
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(a) Optical flux as function of 
total mid-infrared emission 
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(b) Optical colours as function of 
total mid-infrared emission 
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Fig. 12. (a): Optical source fluxes (circles - U-band, squares - R-band & triangles - z'-band) as a function of the total mid-infrared 
flux (the total flux in the Spitzer IRAC 3.6, 4.5, 5.8, 8yum, AKARI SI 1, L18W & Spitzer MIPS 24pm bands) for the sources in the 
FU-HYU-GOODS-N field, (b): Optical colours (optical flux/Total MIR flux) as a function of the total mid-infrared flux for the 
sources in the FU-HYU-GOODS-N field. 
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(b) UV colours as function of 
total mid-infrared emission 
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Fig. 13. (a): GALEX UV band source fluxes as a function of the total mid-infrared flux (the total flux in the Spitzer IRAC 3.6, 4.5, 
5.8, 8^m, AKARI SI 1, L18W & Spitzer MIPS 24pm bands) for the sources in the FU-HYU-GOODS-N field, (b): GALEX UV band 
colours (UV flux/Total MIR flux) as a function of the total mid-infrared flux for the sources in the FU-HYU-GOODS-N field. 



the passage of the 9.7 silicate absorption feature through the 
MIPS24 band (see Section|53]l. 



5.5. Silicate Break galaxies in the GOODS-N field 

The FU-HYU colour-colour diagrams can also be used to search 
for interesting populations on the basis of unique colours or 
colour separations. As an example, we use the colour infor- 
mation in our catalogue to segregate high absorption dusty 
galaxies from the general population. Takagi & Pearson (2005) 
postulated that dusty, heavily obscured starbursts could be de- 
tected via their deep silicate absorption feature at 9.7pm us- 
ing a "drop-out" technique. In this technique, a sharp dip is 
observed in the mid-infrared colours of these sources as the 



silicate absorption feature is redshifted through one of the 
bands. Takagi & Pearson (2005) coined galaxies selected in this 
way "Silicate Break" galaxies and such galaxies were subse- 
quently identified by Charmandaris et al.(2004) using the IRS 
peak up imager on Spitzer. Takagi & Pearson (2005) performed 
simulations for the detection of potential Silicate Break galax- 
ies for a combination of AKARI & Spitzer filters concluding 
that such sources would exhibit colour drops due to a deficit 
in their MIPS 24pm flux between redshifts of 1.2 - 1.6. In 
Figure \T8\a) the MIPS24/L18W colour is plotted against the 
L18W/IRAC8 colour with the galaxy template spectral en- 
ergy distributions for a moderate starburst and ULIRG over- 
laid. The passage of the silicate absorption feature through 
the mid-infrared bands can be clearly seen in the starburst 
and ULIRG templates and map a distinct tail towards lower 
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Fig. 14. Colour-Colour diagrams for sources in the FU-HYU-GOODS-N field as a function of the total mid-infrared flux (the total 
flux in the Spitzer IRAC 3.6, 4.5, 5.8, 8/zm, AKARISll, L18W & Spitzer MIPS 24/mi bands), (a): as a function of the Spitzer IRAC 
8yum band, (b): as a function of the AKARI IRC S 1 1 band. 
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Fig. 15. Colour-Colour diagrams for sources in the FU-HYU-GOODS-N field as a function of the total mid-infrared flux (the total 
flux in the Spitzer IRAC 3.6, 4.5, 5.8, 8/mi, AKARISll, L18W & Spitzer MIPS 24pm bands), (a): as a function of the AKARI IRC 
L18W band, (b): as a function of the Spitzer MIPS 24pm band. 



MIPS24/L18W colours from -0.15> lg(MIPS24/L18W)>-0.4 in 
the colour-colour plane. In Figure [T87a). we adapt the simu- 
lations of|Takagi & Pearson (2005)] for the AKARI L18W band 



plotting the MIPS 24/L18W colour as a function of redshift for 
an ensemble of template star-forming galaxies from the libraries 
of [Tkkagi et al.(2003) (also included are comparative quiescent 
galaxy spectral templates from |Dale et al. (200 1)[ ). Enforcing a 
selection criteria of lg(MIPS 24/L18W) colours < -0.15 should 
preferentially select galaxies in the redshift range of 1.2 - 1.6. 
In FigurefTSTc) the lg(MIPS24/L18W) colour is plotted against 
the MIPS 24pm band flux for the FU-HYU sources. The thresh- 
old for the Silicate Break detection is marked as a solid line 
and it can be seen that there are approximately 9 potential 
candidates in our sample with faint MIPS 24pm fluxes and 



lg(MIPS24/L18W)<-0.15. Extracting all available photometry 
for the 9 candidate Silicate Break galaxies allows the creation of 
photometric spectral energy distributions shown in FigurefTSTt/). 
The presence of the strong emission (1.6/mi bump and 7.7pm 
PAH) and absorption (9.7yum silicate) features can be clearly 
seen in the SEDs. We have used the photometric redshift code 
of Negrell o"et al. (2009)| to fit model spectral energy distribu- 
tions to the 9 Silicate-Break galaxy candidates in Figue[l9] The 
model parameters used for the fits, and the resulting redshift and 
luminosity contributed by the starburst and AGN components 
respectively are tabulated in Table [5] where all the quoted errors 
represent the 99 per cent confidence interval. In total 8/9 of our 
candidates are well fitted by SEDs with silicate absorption in 
the redshift range 1.2 - 1.6, i.e. the expected range sensitive to 
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the Silicate-Break drop out method. The erroneous candidates, 
ID933 appears to have a photometric redshift of z=0.84 assigned 
although it seems plausible that even in this case the 1 .1pm 
PAH feature is being mistaken for the 1 1 .2pm PAH feature and 
that in fact this galaxy too could lie in the preferred redshift 
range.. The brightest source from Figure [19] (ID 1772) is in fact 
the ISO-HDF source HD F PM3 3 (ISOHDF3 J123634+621238, 
|Goldscmidt et al.(T9 97)) at a redshift ~1.28. The galaxy is a 
dusty ultra-luminous merging sub-millimetre source with strong 
emission features and a strong silicate absorption feature in its 
mid-infrared spectrum (Fray er et al.(2008)] [Pope et al. (20 08)). 
The photometric redshift is correct to within 10%, highlight- 
ing the successful identification of dusty galaxies in the silicate- 
break redshift range. 

6. Summary 

We have presented the initial results from the AKARI FU-HYU 
Mission Program in the GOODS-N field describing the data 
reduction process including the additional processing steps re- 
quired to analyze the FU-HYU data due to the omission of 
dithering cycles when the observations were originally taken. 
Combining all the available data in the GOODS-N field, a fi- 
nal FU-HYU catalogue has been produced containing more than 
4393 sources with almost 200 sources detected in the AKARI 
bands. Using the combination of the AKARI and Spitzer multi- 
wavelength coverage a total mid-infrared flux has been defined 
to be representative of the total mid-infrared emission of the 
sources and a measure of the total luminosity of the galax- 
ies (c.f. |Elbazetal. (2002)| >. This total MIR flux is tightly cor- 
related with the longer wavelength mid-infrared (IRC L18W, 
MIPS 24) bands but less correlated with the shorter IRAC 
and IRC bands. The implication is that the Spitzer MIPS24 or 
AKARI L18W band are representative of the total mid-infrared 
luminosity of these galaxies in agreement with the results of 
|Bavouzet et al. (2008) who have indicated that the 24pm flux is 
a good indicator of the mid and bolometric infrared luminosity 
of star forming galaxies. 

The mid-infrared colours have been used to track the pas- 
sage of the PAH emission features through the observation bands 
and to segregate the star-forming population at z~ 1 which have 
higher total MIR / single band colours than the more quiescent 
population. In particular we have shown that an excess in emis- 
sion in the AKARI S 1 1 band is indicative of a moderate redshift 
population. 

Using the AKARI IRC LI 8 W to Spitzer MIPS 24 band colour 
as an example diagnostic, we have shown that it is possible to 
segregate specific populations in the colour colour plane, thus 
we have used the "Silicate-Break" technique to extract extinct, 
dusty galaxies from our FU-HYU sample. This population is 
sensitive to the passage of the 9.1pm absorption feature through 
the MIPS 24 band in the redshift range -1.2-1.6 and this tech- 
nique has been successful in identifying 8 possible candidates in 
the GOODS-N field. Models fits of spectra using a photometric 
redshift code have indeed confirmed that the sources lie in the 
redshift range expected for Silicate-Break galaxies confirming 
that the silicate break selection method can provide a powerful 
means to detect dusty ULIRGs at moderate redshift. 
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Fig. 17. Infrared colour-colour diagrams for sources in the FU-HYU-GOODS-N field as a function of the AKARI L18W (or SI 1) 
band and the total mid-infrared flux (the total flux in the Spitzer IRAC 3.6, 4.5, 5.8, 8jum, AKARI SI 1, L18W & Spitzer MIPS 
24pm bands), (a): as a function of the Spitzer IRAC 8/mi band, (b): as a function of the AKARI IRC SI 1 band, (c): as a function 
of theSpitzer IRAC 8pm band (to AKARI IRC S 1 1 band), (d): as a function of the Spitzer MIPS 24pm band. Model predictions for 
colour-colour tracks are also shown for normal galaxy (dashed), star-forming galaxy (dotted), ULIRG (triple dot-dash) and AGN 
(solid) spectral templates. The large squares mark the zero redshift template and the crosses along the tracks are spaced at equal 
redshift intervals of Az=0.2. 
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Fig. 18. Selection of potential Silicate Break Galaxies in the FU-HYU sample, (a) The MIPS24/L18W plotted against the 
L18W/IRAC8 colour. Overlaid are the SEDs of a moderate starburst and an ultra-luminous galaxy. The SED tracks are marked 
by a large square at the redshift =0 position and have crosses a step of 0.2 in redshift thereafter extending to a redshift of 2. The 
passage of the silicate absorption feature through the mid-infrared bands can be clearly seen in the starburst and ULIRG templates 
and map a distinct tail towards lowerMIPS24/L18W colours in the colour-colour plane, (b) MIPS24/L18W mid-infrared colour as a 
function of redshift for a variety of star-forming galaxies adapted from Takagi ^fe Pearson (2005)| Dotted lines ar e the starburst mod- 
els of Takagi et al.(2003) solid lines are comparative quiescent galaxy spectral templates from |Dale et al. (200 1)[ The silicate break 
selection criteria is set to be ~lg(MIPS24/L18W)<-0. 15. (c) The diagnostic lg(MIPS24/L18W) colour against the MIPS 24/zm band 
flux. The solid horizontal line represents the Silicate Break detection criteria of lg(MIPS24/L18W)<-0. 15 below which around nine 
potential candidates are selected, (d) Mid-infrared photometry in all available bands for the nine candidate Silicate-Break galaxies 
in the FU-HYU GOODS-N field. Several emission features are highlighted in the observed frame of the galaxies implying that all 
(except one) sources are between redshifts 1.2 - 1.6. 
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Table 5. Best fit results for the 7 objetcs in Figure[19]for candidate Silicate Break galaxies. 
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Fig. 19. Model fits for the candidate Silicate Break galaxies using the photo-metric redshift code of Negrello et al. (2009) 



The dashed and dotted lines are the fitted starburst and AGN components from the spectral libraries of Takagi et al.(2003) & 
Efsta thiou & Rowan-Robinson (19 95) respectively and the solid line is the composite fit. Photometric redshift fitting parameters are 
listed in each panel and tabulated in Table [5] following Negrello et al. (2009) 



